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ABSTRACT 


High resolution (10 km s“^) [N II] echelle spectra, sampled every 13 arcsec 
in a 6' X 6' region around SN1987A were obtained on the CTIO 4m telescope. 
The map shows a complicated velocity structure consistent with that reported 
previously for the interstellar medium (ISM). Three components, Vhei = 265, 277 
and 285 km s“^ were identified as N157C (or the R1170-complex: |Xu, Crotts| 
fc Kunkel 1995|) . The radius of this large superbubble was found to expand at 
10 km s“^, with a lifetime of 6 x 10® years and a total energy of 3 x 10®^ ergs 
determined from its radius and velocity according to superbubble theory 
( [McCray fc Kafatos 1987] ). The Vhei = 235 km s“^ component corresponds to 
the near side of 600 pc giant superbubble reported earlier. This bubble is over 
10^ years old, and has blown out of the LMC disk. Two other components, 
Vhei = 255 and 245 km s“^ are identified as the inner major light echo ring (a 
double-shell structure) at about 130 pc in front of SN1987A. There are also two 
high velocity components, 300 and 313 km s“^, which are possibly the far side of 
a superbubble in which SN1987A exploded. We also noticed that there are two 
components at 269 and 301 km s“^ within 20" of SN1987A. These structures are 
probably due to the emission from the progenitor star’s red supergiant wind. 
We find that the time it took the SN1987A progenitor to move to the current 
location 300 pc behind N157C is comparable to the age of N157C as well as that 
of the progenitor itself. 


Subject headings: galaxies: Magellanic Clouds - interstellar: matter - stars: 


super novae 
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1. Introduction 


30 Doradus in the Large Magellanic Cloud (LMC) exhibits extremely active star 
formation and interstellar activity. It serves as an interesting laboratory for the study of 
the evolution of the interstellar medium (ISM) and its interaction with massive stars. Many 
giant gas sheets in this region were identihed in deep Ha photographs ( [Meaburn 1980|) 
and in 21-cm emission (jMeaburn et ah 1987|) . These shells enclose about ten interlocking 
superbubbles each about 100 pc in diameters and expanding at about 30 - 50 km s“^. 
Emission and absorption lines ( [Smith fc Weedman 1971| ; [Canto et al. 1980| ; [Cox &] 
Deharveng 1983 ; [Meaburn 1984] ; and Clayton 1987[) revealed a complicated velocity 
structure which closely correlates with HI components. The observed components range 
over Vhei = 185 - 307 km s“^. X-ray maps ([Wang &: Helfand 1991a[ , b) demonstrate a 
good correlation with Ha images, indicating that these superbubbles contain hot gas. It 
was shown ([McCray fc Kafatos 1987[) that repeated supernova explosions of massive stars 
play the major role in pressurizing these giant bubbles and evacuating interstellar space, 
culminating in 600-800 pc diameter supershells. 

SN1987A in the LMC provided a unique opportunity for studying part of the 
30 Doradus complex. Ca II, Na D and UV absorption lines against SN1987A ( Andreanp 


Eerlet fc Vidal-lVIadjar 19871 , [Vidal- Aladjar et al. 1987[ , [Alagain 1987| , and [Blades et al. 1988]) 
exhibit groups of LMC components near Vhei = 129, 171, 206, 226, 250-265 and 281 km s“^. 
Other studies ([deBoer et al. 19^ ) indicate that the components with velocities less than 
140 km s“^ are not associated with the LMC (and are also weaker). 


The study of this part of space was greatly boosted by the discovery of the light echoes 


from SN1987A ( [Crotts 1988[) . As forward scattered light from the SN, light echoes carry 
information not only about the SN, but also about the interstellar dust which reflect them. 
In the decade since the hrst discovery, these light echoes from SN1987A have revealed 
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a great deal of information about this region ([Suntzeff et al. 1988 ; Gouiffes et al. 198S ; 
[Crotts, Kunkel fc McCarthy 1989|; Pond et al. 1989| ; [Couch, Allen fc Malin 1990| ; pCu, 


Crotts fc Kunkel 1994|; and ^pyromilio et ah 1995]) . In a more recent paper (jXu, Crotts fc 


[Kunkel 1995|) , we discuss comprehensively the large scale ISM structures in the SN1987A 
foreground, and constructed a three-dimensional (3D) map of these complicated structures 
from the light echoes. Twelve dust sheets were identified: R310, R430, W700, S730, N980, 
R1170, SW1400, NW1500, W1700, R1830, SE3140, and N3240. These names are derived 
from the following rules: the hrst one or two letters indicate the direction relative to 
SN1987A, e.g. “R” for ring shaped, “N” for north and “SE” for south-east; the following 
number is the distance in front of SN1987A in light years. These structures enclose a 110 pc 
radius superbubble which interlocks with a smaller, 30 pc bubble, a 600 pc diameter giant 
superbubble, and a shell which likely envelopes SN1987A. The giant superbubble appears 
to have also been observed by another group ([Bruhweiler, Fitzurka fc Gull 199l[) as an 
oval-shaped dust ring of 1.5-2.0 kpc in dimensions in a narrow continuum band (4215A). 

One of the important results from 3D light echo mapping is that the superbubble 
N157C is not just the smallest “cavity” enclosed by the Ha hlaments to the north-east 
of SN1987A (for example, see [IVfeaburn 1980[) . Instead, the complicated Ha structure 
which extends to the south-west of SN1987A, has a spherical 3D structure whose center 
of curvature falls on LH90. Therefore, we think that the N157C supershell extends to a 
radius of about 115 pc to the south-west, in the direction of SN1987A. An X-ray map 
( [Wang &: Helfand 199 la[) reveals a bubble bigger than the small cavity which was believed 
to be N157C. The X-ray maps are likely to show where gas is hitting the inner surface of 
the shell, with gas interior to this probably too diffuse to detect. The lowest contour level 
(0.0013 photon s“^ arcmin”^) of the partial X-ray ring around the central OB association 
LH 90 actually extends to a radius of about 60 pc (4.0 arcmin), but is incomplete to the 
south and west around its circumference. It appears that the N157C bubble extends to 
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larger radii towards the west and south, towards the direction of the SN. In the deepest 
ROSAT PSPC images of the same held (pdasinger et ah 1996|) , it is clear that two segments 
of X-ray emitting gas, roughly concentric with respect to the center of N157C, sit in 
the region extending into our survey held, the farthest one at a radius of about 100 pc 
(~ 6.5 arcmin) with respect to N157C’s center. Future observations might reveal if these 
have the same spectrum as the rest of the X-ray shell. One cannot prove on the basis of 
X-ray morphology alone where the edge of the superbubble extends in the direction of the 
SN, but it seems likely to be larger than the 60 pc radius seen at other position angles. In 
contrast N157C seen in 21-cm emission from H I extends over at least 200 pc in diameter 
( [Kim et ah 1998|) , nearly coextensive with the dust echo. 

It is highly desirable to study the velocities of these structures. Echelle emission 
spectra in the vicinity of SN1987A were reported before ( jMeaburn 1990|) . The components 
were found consistent with absorption structures, with 255 and 280 km s“^ the two 
most ubiquitous components. In this paper, we will present the results of our large held 
echelle survey around SN1987A and its relation to the ISM structures. Section 2 discusses 
observations and data reduction. Section 3 will explain our numerical model and htting 
procedure. We interpret the velocity map in §4. 


2. Observations and Data Reduction 


The survey was conducted on 22-24 January 1993 and on 16-18 January 1994 using the 
Cassegrain echelle spectrograph on the CTIO 4m telescope. We used an echelle setup of hve 
evenly spaced long slits, oriented north-south. The slits were 131" long and centered 12".5 
apart from each other. With such hve slits, every exposure covered a held of 131" x 50". 
By making the ohset between helds 120" north-south and 60" east-west (hence 10" between 
adjacent exposures) there was sufficient overlap north-south to check that we had returned 
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to the same right ascension on later nights (while the east-west offsets, also donble-checked 
by occasionally re-acqniring the SN, were made in qnick snccession on the same night). 
We planned to nse 18 snch evenly distribnted fields to cover a 370" x 350" region, from 
125" sonth to 245" north of SN1987A, and from 170" west of to 180" east of SN1987A. We 
managed to get 17 fields, not having enongh time to observe the sonth-westernmost corner 
field. Additionally, two single 247"-long slits oriented east-west, one centered on the SN and 
one centered 100" north of the SN, were also nsed beforehand to check for strnctnre that 
might extend beyond the 120 km s“^ free velocity range between adjacent slits. No snch 
strnctnre was fonnd, with the exception of a small region in the five-slit data at the edge 
of the field sontheast from the SN. This corresponds to the “Honeycomb Nebnla” ( |Wang] 
19921 , pVleabnrn et al. 1993| , [Meabnrn et ah 1995| , |Chn et ah 1995|) , where we ignored some 


data which appeared snsceptible to overlap. The confignration of all of these slit positions 
is shown in Fignre 1. 


We chose to look at the [N II] donblet (air A6548.08 and A6583.41: Psterbrock, Tran, 


fc Viellenx 1992|) , and nsed an Ha filter (central wavelength = 6563A, FWHM=75A) to 
separate orders. We nse [N H] instead of Ha becanse the thermal width of Ha line is abont 
3.7 times of that of [N H] line at the same temperatnre dne to the difference in their atomic 
mass, a critical difference. Resolntion of the setnp is 10 km s“^ near [N H] A6583.41, and 
the free velocity range between every two slits is 120 km s“^. Exposure time of each field is 
3600 seconds. 


We used IRAF to reduce the data. The process started with column bias, two- 
dimensional bias frame subtraction, and dark current subtraction. Bad pixels were 
interpolated over. The response function image from each night, which was later used to 
fiat field the images, was constructed from combined quartz exposures using the IRAF 
“response” task. We then used the IRAF “illumination” task to fit “slit functions” from 
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sky exposures. These slit functions were then applied in illumination correction to all the 
object images to give a flat response along slits. We calibrated the spectra into wavelength 
space using ThAr comparison spectra. The calibration error can be assessed by comparing 
the measured sky line wavelength to the standard values ( [Hubbard et al. 1986| ). We found 
the error was less than 0.04 A, or equivalently 2 km s“^ at A6548. This error is only 
1/5 of the 10 km s“^ resolution. Because the later process requires a high S/N ratio, 
we binned the spectra over every 13".3. A continuum was then £t and subtracted from 
each one of the 13".3 spectra. We then scaled the spectra in adjacent helds to the same 
flux level. The scaling was possible because we deliberately overlapped each held with its 
neighbors. We multiplied spectra of each held by a constant so that neighboring helds have 
the same Ha and [N II] (A6583) hux in the common region. Finally, night sky lines were 
properly subtracted from the spectra. Unfortunately, one of the [N II] doublet (A6548.08) 
was terribly contaminated by OH-lines. Therefore, we only look at the other [N II] line 
A6583.41. 


3. Velocity Components Analysis 


We transformed the spectra into velocity space using the Doppler shift formula. The 
[N II] vacuum wavelength are A6549.87 and A6585.21 which are calculated from their air 
wavelengths with Edlen’s formula (jEdlen 1966|) . 


3.1. Component Identification 

It is crucial to first identify structures corresponding to each slit, because the 5-slit 
echelle may cause cross-slit overlap. For this purpose, we looked at two single slit spectra: 
one centered on SN1987A, and the other about 100" to the north. We hrst binned them 







over 130" to achieve a high S/N level. Figure 2 show Ha and [N II]A6583 lines of these 
spectra. In Panels 1 and 3 (Ha line), the two structures near 0 and 650 km s“^ are skylines 
A6562.67 and A6577.21 ([Hubbard et al. 19861) respectively. In Panels 2 and 4 ([N H] A6583 
line), the structure near 600 km s“^ is a skyline A6596.55 ([Hubbard et al. 1986|) . These 
spectra show that the only interesting structures are in 210-330 km s“^ interval, which is 
within our slit-to-slit free range. Although this conclusion is technically correct only in the 
region on and to the north of SN1987A, we think it is very unlikely that some extraordinary 
low or high velocity structure would appear only to the south of SN1987A, since most 
structures are seen to be weaker there. 


3.2. Numerical Analysis 

The goal of numerical analysis is to identify the number of velocity components and to 
measure the properties of each component, e.g. velocity, intensity and FWHM. Our analysis 
was a two-step process. We hrst £t the data with various numerical models with up to five 
gaussian components. Then based on the htting result, we used a statistical algorithm to 
select the correct component numbers. 

First, we varied the number of Gaussian components from one to hve to produce hve 
different models. These models were htted separately to the spectra using least htting to 
determine the parameters. Because the temperature of each component in the LMC should 
not change very much, and hence the FWHM of each component should be roughly the 
same, we forced a uniform FWHM across components within each model, while the height 
and position were allowed to vary freely and independently from each other. To check out 
the “uniform FWHM” assumption, we repeated the whole process but let FWHMs vary 
as free variables. The “model probability” method (see below) selected the same number 
of components in every held, regardless the assumption on FWHMs. We also noticed that 
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the velocities of corresponding components are a little different, but within our resolution. 
Therefore, we think our modeling based on the “uniform FWHM” assumption is robust. 
The htting was performed with the “deblend” option of IRAF’s “splot” task. We chose to 
work interactively because the components are normally not well separated and automatic 
htting tends to yield wild results. 

Based on the least htting result, we computed the probability that a “correct” 
model would yield a higher than the actually htted result. The selection rule is that the 
most likely model has the maximum probability. It has been shown that the probability 
dehned in this way satishes a x^-distribution and consequently can be calculated ( |Fres^ 
et al. 199^) . For every prohle, we compared the hve models based on their probabilities, 
and select the one with the highest probability. In most cases, we found the “correct” 
model yields a probability several orders of magnitude higher than other models, and 
hence convincingly establishes itself. For a few prohles, two or three models have similar 
probabilities. In these situations, we made our choices by continuity argument. We look at 
the neighboring locations for any dehnitely “correct” models, and select the most similar 
one. Figure 3 shows the reduced data from one of our 17 helds, detailing the [N II]A6583 
line, and the result from our velocity component analysis graphed on a corresponding scale 
in the adjacent diagram. The analysis reveals all of the major features and nearly all weaker 
ones, and follows gradients in the velocities of components with reasonable accuracy. 

With the aforementioned process, we determined at every grid point over all 17 helds 
the number of components and properties of each component. We then looked into the 
connection among these grid points for collective structures across the held. The result will 
be discussed in the next section. 
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4. ISM in the Foreground of SN1987A 


This survey shows the following ionized clouds: Vhei = 237, 245, 255, 265, 277, 285, 300, 
and 313 km s“^. (When clouds are within about 5 km s“^ of each other, they tend to blend 
together.) These are shown in Figure 4a. The 265 — 285 km s“^ grouping is sufficiently 
complex that we show it separately in Figure 4b. Four of these components (255, 277, 300 
and 313) were also reported in a smaller field containing SN1987A (jMeaburn 1990|) . Two 
of them, Vhei = 255 and 277 km s“^ were discovered in the interstellar absorption spectra 
toward SN1987A (see the citations in §1). A HI 21-cm study in this region exhibited similar 
profiles: Vhei = 243, 274, 296 and 318 km s“^ ( [Meaburn, McGee fc Newton 1984| ). Similar 
structures were also discovered in the halo of 30 Doradus from [O III] lines ([Meaburn 
|1984| ; [Meaburn 1988]) . In our survey, we found the 245 and 255 km s“^ components 
(Group B in Figure 4a) have similar properties and hence are likely to be two layers of one 
structure (see §4.3 below). Similarly, Vhei = 265, 277 and 285 km s“^ clouds (Group G and 
Figure 4b) are seen to be interrelated with each other and even connected in some places. 
(Some components are seen to disappear into others when approach within the 5 km s“^ 
blending range.) We think these three are really one structure (see §4.1 below). The two 
high velocity clouds (300 and 313 km s“^; Group D) are likely to be the far side of the 
245-255 km s“^ shell (§4.3, see also [Meaburn 1990|) . The other component (237 km s“^; 
Group A) is reported for the first time in this study. It is interesting that we did not see 
this component in the two single long slit spectra, one centered on SN1987A and the other 
100" north of the supernova. We show in §4.2 that this component correlates to the two 
clouds about 1 kpc in front of SN1987A (PCu, Grotts fc Kunkel 199^ ). 

















4.1. N157C 


The Vhei = 265, 277 and 285 km group (C in Figure 4a, also seen in Figure 4b) 
is the most far reaching structure in this held. They not only are interrelated with each 
other, but also have comparable hux level consistently bright throughout the whole held. 
Similar behavior was observed in [O III] and Ha emission lines (jMeaburn 1990|) , in the HI 
21-cm survey ([Meaburn, McGee fc Newton 198'^) , and in Ca II, Na D and UV absorption 
lines (see §1 for references). Therefore, it is tempting to recognize these three components 
as one dominant structure in front of SN1987A. In an earlier paper (pCu, (Jrotts fc Kunkel 


1995|), we show that the most prominent structure in the SN foreground is the echo 1170 


light-years in the SN foreground, the R1170-complex (identihed as N157C). This similarity 
marks the hrst correlation indicating that these velocity components are indeed N157C. 


As another piece of evidence, we plotted the hux contour of these components on top 
of the 2D image of the ISM dust structure, the R1170-complex ( |Xu, Crotts fc Kunkel 19951) 
(Figure 5). This plot shows a general tendency that both structures are bright to the north 
of SN1987A. In addition, we observed a good correlation not only on the large scale but 
also at many small bright spots, limited by the spatial accuracy of the echelle survey of 13" 
east-west and 40" north-south. 


We also observed that the two structures split in the same way. The 3D light echo 
map (PCu, Crotts fc Kunkel 19951) demonstrated that the cloud, consisting of N1500 and 
W1700, splits from R1170 to the north of SN1987A and joins R1830 southwest of the SN. 
Similarly, Figure 4b shows the 277 km s“^ component joining the 285 km s“^ component to 
the south (and slightly east) of SN1987A and the 265 km s“^ component to the southwest. 
This geometry is consistent with the picture that the 285 and the 265 km s“^ components 
are the far and near shell of N157C, while the 277 km s“^ component is a cloud joining the 
two structures to form a superbubble. 
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N157C, with a diameter of about 200 pc, is expanding at a speed of 10 km s“^, 
therefore. We believe this bubble is still in the pressure-driven stage for the following 
reasons. First, this size is common in other parts of the LMC, particularly in 30 Doradus 
( [Meaburn 1980|) . Most of these giant bubbles are young and rapidly expanding ([Meaburn et 
|al. 19871) . McCray and Kafatos (1987) attributed the generally large superbubble size to the 
thick gas disk of the LMC and its low metallicity. Secondly, the interior of a pressure-driven 
bubble is hlled with hot gas as opposed to the cool interior in the snow-plow stage. In 
addition, Wang and Helfand (1991b) reported X-ray (~ 0.8 — 1.2keV) emission in N157C, 
corresponding to very hot gas in its interior. Therefore, N157C is in the energy-conservation 
stage. McCray and Kafatos (1987) gave the radius and velocity properties for these bubbles: 


Rs = 97 pc 

(1) 

Vs = 5.7 km s-^ 

(2) 

Dividing one equation by the other. 


^ Rs 5.7 km s~^ 

Vs 97 pc 

(3) 


where uq is the gas atomic density, and tj = t/l(f yr. Let Rs = 100 pc, 1^ = 10 km s“^. 
We estimate the age of N157C to be f = 6 Myr. Therefore, this superbubble is still very 
young. It is estimated ([McCray fc Kafatos 1987D that stellar winds from a rich cluster 
(iV* > 10) dominate its superbubble expansion for a few million years until supernova 
explosions become more important, while the superbubble may expand to less than ps 100 
pc. We think N157C has just passed the stellar-wind-driven stage, and is in rapid expansion 
driven by powerful supernovae from LH90. 

From other sources, we estimate the neutral hydrogen column density Nhi of the 


supershell N157C. From 21 cm emission ([Luks fc Rohlfs 1992|) , about 70% of the H 1 
column density towards LH90 is contained in a narrow component centered at 270 km s“^ 
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with a total column density of all components of about 2.2 x 10^^ cm“^. This total Nm is 
corroborated from the pointing (^^43) containing N157C from the survey of |Meaburn et| 
|ah 1987| : Nhi = 2.3 x 10^^ cm“^, and the value inferred from the damping wings of the Ly 
a absorption line towards the star Sk —69°203, in the direction of N157C: 2.4 x 10 ^^ cm ^ 
( [Fitzpatrick fc Walborn 19^) . Thus the value of N^i for N157C is about 1.5 x 10^^ cm“^. 
(However, none of these surveys have resolution better than about 15', so we cannot be sure 
of the effects of substructure.) Since this column density is optically thick at the Lyman 
limit, we will assume the gas is shielded and therefore largely neutral. Accounting for 
projection effects (a factor of four), the kinetic energy of the expanding HI shell is therefore 
Esheii = 5.4 X 10^° ergs. This implies a total kinetic output from SNe that is 77/15 times 
greater ([Mac Low &: McCray 19881) : Es = 2.7 x 10®^ ergs. Each typical Type H supernova 
of a massive O star releases about Esn = 10®^ ergs of kinetic energy. Given an average 
supernova rate of about one per million years in a typical OB association ([Mac Low & 
McCray 1988[) , the energy budget is consistent with the aforementioned age estimate. 


Given Egheii determined from theoretical consideration of a SN driven expansion 
( [McCray fc Kafatos 1987[) : 


Esheii = 4.0 X lO^^er^is (^'*^ 51)^7 


(4) 


we hnd ALE 51 = 22. Combining this with Equation 1 (or 2) yields an initial density no = 4. 
This is a reasonable estimate, corresponding to the Nm = 2.4 x 10^^ cm“^ spread over a 
scale height of 100 pc, which is reasonable in terms of the observed 3-D structure. 

On the other hand, from stellar content point of view, the aforementioned estimate 
Ni,E^i implies several massive SN explosions in LH90, given that the typical energy of a 
Type H SN, Esi, is a few for the most massive stars (M > 30Mq). The main-sequence 
lifetimes of massive stars are given approximately by tus ~ 30 Myr (M*/[10Mq])““, 
where a 1.6 for 7 < M* < 30M* ( [Stothers 1972[ ) and by fjus ~ 9 Myr (M*/[10Mq]“°-^ 
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for 30 < M* < SOM* (|Chiosi, Nasi fc Sreenivasan 1978|) . In a time scale of 6 Myr, only 
a few very massive stars with (M*>30Mq) would explode. It is estimated that an OB 
association should produce roughly 9 times as many stars with masses in the range 7-30 Mq 
(main-sequence spectral type B3-08, corresponding to iV*) as stars with masses > 30Mq 
(MS type 07-03) ( [McCray fc Kafatos 1987|) . Given that there are about 50 stars with 
mass > 7Mq ([Hill et ah 19941) , there should be about 5 massive stars of M*>30Mq which 
should have exploded. 


We can also compare our predicted SN explosion number to empirical determinations 
of the stellar population. In addition to the original study of LH 90 (p.jucke fc Hodge 1970|) , 
two recent studies of the stellar population of the association yield somewhat different 
results. Recent UIT data ([Hill et ah 1994|) on LH90 and six other LH OB associations in 
all but the inner portion of the 30 Dor complex give the number of stars and their inferred 
mass. For the clusters with the largest mass stars (~ 50Mq), the ratio of M > 30Mq versus 
7Mq > M < 30Mq stars predicted in the previous paragraph from theoretical estimations 
appears inconsistent with the UIT data, in the sense that more massive stars are actually 
seen by UIT. From these data, LH90 appears to have stars as massive as 41Mq, while LH89 
and LH99 show about eight stars more massive than this (out of a total of 40 stars in LH90 
and 68 stars in LH89 and LH99 more massive than 15Mq). Thus LH90 might be missing 
approximately hve stars of M Ri 50Mq (or more if stars have exploded in LH89 or LH99); 
these have presumably formed SNe. 


We should note, however, the existence of stars in LH90 more massive than those 
allowed (> 30Mq) by a 6 Myr year turnoff age. Perhaps some massive stars were created 
millions of years after the hrst stars. Indeed, our prediction for the number of massive stars 
seems more consistent with the results from the other study ( [Festor et al. 1993[ ), which 
uses optical photometry and spectroscopy of LH 90 to derive spectral types and Mf,oi values 
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consistent with abont 7 stars with M* >30Mq versus more than 47 stars with mass > 7Mq 
(according to our simple conversion of Mi,oi to M*). Indeed, their study concludes that the 
association contains a sub-population 3-4 Myr old and a possible larger sub-population 
7-8 Myr old (P?estor et al. 199^) . This is consistent with our arguments regarding the SN 
rate. 


4.2. 600 pc Giant Superbubble 


Xu, Crotts and Kunkel (1994) hrst reported two dust bodies at about 1 kpc in front of 
SN1987A. They (pCu, Crotts fc Kunkel 199'^) also found that the southeast piece of these 
two correlates with one of the bright Ha filaments extending from N157C discovered by 
Meaburn (1984). They went one step further to suggest this filament as the near side and 
N157C as the far side of a giant superbubble of 600 pc in diameter. Consistent with this 
picture, a component at 235 km s“^ (“A” in Figure 4a) has been found in the regions 
far away from SN1987A. Figure 6 plots flux contour of this component over the image of 
SE3140 and N3240 (pCu, Crotts fc Kunkel 19951) . We found that one piece of the 235 km s“^ 
components corresponds exactly to N3240, and another piece correlates with SE3140. It is 
unfortunate that we only had enough time to finish the survey to the north of Dec —69° 18' 
and to the west of RA 5h 36m 00s (in J2000) delineated by the dashed lines. Still, the 
available data indicates that 235 km s“^ gives a good correlation with SE3140, expanding 
at about 40 km s“^ relative to the centroid of N157C. 

It is interesting to see the 235 km s“^ component distributed in a patchy ring-like 
shape. Recent light echo images (1995-1997) showed N3140 becoming larger and extending 
to the position angle PA = 45° where the (5h36m, —69° 14') 235 km s“^ components is 
located. Therefore, we think the 1 kpc ISM structure may indeed show a hole in front 
of SN1987A. Moreover, we have several reasons to believe this hole may relate to the 
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blow-out of the 600-pc-diameter superbubble. First, we have not observed any structure 
within 3' of SN1987A either in the light echoes at 1 kpc or in the echelle survey at 
235 km s“^, indicating that the cloud has been fragmented. This result can be explained 
by a superbubble blow-out. Secondly, Wang and Helfand (1991b, c) reported that this 
region has soft X-ray emission much higher than the background but much lower than other 
young superbubbles in 30 Doradus. This region may even be contained in the possible 
giant superbubble LMC-2, roughly 1 kpc across ( |Wang fc Helfand 1991c|) which has blown 
out of the disk. Finally, 600 pc is one of the largest superbubbles ever observed in the 
LMC (pVleaburn 1980|) , and therefore may represent a size close to the upper limit imposed 
on the superbubble size by the LMC disk thickness. In conclusion, we think this 600 pc 
diameter giant superbubble appears to have blown out of the LMC disk and is expanding 
in radius (assuming it is pushing on one side against N157C) at 20 km s“^. McCray and 
Kafatos (1987) suggested that when the superbubble cools or breaks out of the galaxy, it 
will be distorted by gravitational instability, collapse and form giant clouds, and hnally 
ignite propagating star formation. LH90, which formed N157C, may have itself formed in 
this way in the gravitational collapse of the 600 pc giant bubble. 

If this superbubble blew out in its snow-plow stage, prior to this it was described in 
terms of the time (tc) and expansion radius {Rc) beyond which cooling becomes important 
in the interior e.g. post-adiabatically according to ( [McCray &: Kafatos 1987]) : 

R{t) ^ Rc{t/Ry^^ (5) 

for 

R, = 50 pc C-°-®(W^5i)°-" no (6) 

and 

4 = 4 X 10 ® C-'-'^(W^5i)°-' 


( 7 ) 
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where ( is metallicity, iV* is number of stars > 7Mq, E^i = E/lO^^ergs, hq is gas atomic 
density. In the LMC, ( ~ 0.3 (|Dufour 1984|) and a density typical for the LMC (several 
times lower than near N157C), riQ ^ 0.35 ( [Hindman 1967| ). Let’s guess Ah = 50 (roughly 
that of N157C) and E^i = 1, hence tc = 160 Myr, and Rc ~ 1.3 kpc, much bigger than 
300 pc. Therefore, we think this giant bubble was probably hot when it blew out, implying 
an age at the time of blowout of about 10 Myr from eq. 3, or 15 Myr using eq. 1 and the (less 
certain) parameters above. These are probably consistent with the age of LMC-2 ( |Wang| 
jfc Helfand 199 Icj) . This appears older than N157C, indicating the 600 pc superbubble 
predates it, and may have even brought about LHOO’s formation. 


4.3. Interstellar Environment of SN1987A 


Xu, Crotts & Kunkel 1995 demonstrated that the two dense clouds (R430-complex 
including R310) about 100 pc in front of SN1987A do not correspond to any prominent 
Ha structures, and therefore are likely related primarily to SN1987A. R430 is brighter 
and larger than R310. The lack of correlation is also consistent with the assertion that 
SN1987A exploded about 300 pc on the far side of the LMC disk (pCu, Crotts fc Kunkel] 
1995]) . Figure 7 plots flux contour of Vhei = 255 and 245 km s“^ (Group B in Figure 4a) 


on top of R430-complex, indicating a good correlation between these two images. On the 
other hand, the 255 km s“^ component is brighter and larger than 245 km s“^ component, 
in good agreement with the relative size and brightness of R430 and R310. The correlation 
between these two studies indicates that the double shell structure R430 and R310 have 
velocities 255 km s“^ and 245 km s“^ respectively. The 255 km s“^ component was reported 
by other groups before (see [Meaburn 1990] and the reference therein). However, Meaburn 
(1990) suggested that the 255 km s“^ shell is within 10 pc of SN1987A. 


Two other ionized clouds Vhei = 300 and 313 km s ^ in the echelle spectra (Group D 
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in Figure 4a) do not correspond to any absorption structures. It was suggested ([Meaburn 
1990|) that these two clouds are likely to be the far side of a superbubble which envelops 
SN1987A. We noticed that this identification is consistent with the double shell structure 
on the near side. Unfortunately, the light echoes have expanded to only several light years 
behind SN1987A. So far, we do not have any direct light echo results on these two possible 
double shells on the far side of the supernova. Therefore, we can only speculate on the 
structures of these two high speed components. 


If these two components are indeed on the far side of a superbubble surrounding 
SN1987A, this superbubble is expanding in radius at 30 km s“^. If we adopt a average 
radius of 90 pc, the age of this bubble is about 2 Myr, (eq. 1). Let the shell density of this 
bubble Nh 1 x 10^° ([Meaburn et ah 1987|) , using the same method in §4.1, we calculate 
the total energy of this superbubble to be 3 x 10^^ ergs, which is the energy of two to three 
typical Type II supernovae. Both the energy budget and the age suggest that there had 
been a small number of supernova explosions before SN1987A in this superbubble. Because 
of the short lifetime of this bubble, these previously exploded supernovae must be very 
massive stars (M* > 30Mq) (jChiosi, Nasi fc Sreenivasan 1978|) . 


The aforementioned superbubble assertion has some problems. First, there are not 
many stars close to SN1987A. Even though one may argue that massive stars might have 
exploded, the lack of low mass stars is still a difficulty. Second, the expansion velocity of 
30 km s“^ is too high compared to that of N157C (10 km s“^) which is pressurized by the 
second richest cluster LII90. The only way around this problem is to assume very low gas 
density in this region. In Equation 2, adopt Vs = 30 km s“^, = 3, = 0.2, we get 

no = 0.02, which is an order of magnitude lower than the typical LMC value of no = 0.35 
( [Hindman 1967]) . Perhaps this is reasonable, given SN1987A’s likely distance from the plane 
of the LMC, but the existence of a superbubble surrounding SN1987A is still not certain. 
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We also noticed that there are two components at 269 and 301 km s“^ within 20" of 
SN1987A. They have much higher flux level than those probably more distant components 
at 245, 255 and 317 km s“^ (and at 300 km s“^ not within 20"). It has been suggested that 
these structures are probably due to the emission from the red supergiant wind (|Crotts| 
fc Kunkel 1991|) . This structure will be treated in a separate paper (jCrotts fc Heathcote| 
19981). These two components have been reported previously ([IVleaburn 1990|). 


5. Discussion 


We note that our light echo 3-D map, in conjunction with the velocities of the matter 
in front of the SN from this study, implies that the SN is well behind most of the material 
in its portion of the LMC, and hints at how this occurred. This refutes claims made on 
the basis of extinction towards a nearby star, Sk —69°203, that the SN is embedded in the 
middle of the main H I complex ( Fitzpatrick fc Walborn 19901) . We approximate the speed 
of LH90 with that of N157C, about 270 km s“^. This should also represent the speed of the 
LMC stellar disk. Therefore, SN1987A is moving away from the LMC disk at a speed of 
20 km s“^ given that the velocity of SN1987A of 289 km s“^ ( |Crotts fc Heathcote 1991| ). 
At this speed, it would have taken the SN1987A progenitor ~ 10^ years to move to its 
current location 300 pc away from the stellar disk. It is interesting to compare this time 
scale to the lifetime of N157C (6 Myr, §4.1) as well as that of the speculated superbubble 
surrounding SN1987A (2 Myr, §4.3), and that of the progenitor itself. We speculate that 
the SN and N157C were created together (followed by the superbubble around SN 1987A), 
and are pursuing further methods to test this. 
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Figure Legends 

Figure 1: A 15' field centered on SN 1987A imaged in a band highlighting Ha and 
[N II], and showing the positions of our echelle long slits relative to features in the field. 
The OB association LH90 sits near the corresponding label, with the most prominent 
groups of stars seen to the west and south of the label. It is surrounded at radii out to 
~ 10' by emission from gas associated with the accompanying supershell N157C. The 
Honeycomb Nebula (labeled “HN”) sits about 2.5' southeast of the SN. At the time of these 
observations, the light echoes extended to roughly the edge of our slit array. 

Figure 2: Binned spectra of two single 120" long slits. From upperleft counterclockwise 
are Panels 1, 2, 3 and 4. Panel 1: Ha spectrum centered on SN1987A; Panel 2: [N H] 
(A6583) spectrum centered on SN1987A; Panel 3: Ha spectrum 100" north of SN1987A; 
Panel 4: [N H] (A6583) spectrum 100" north of SN1987A. The abscissa is Vhei in km s“F 
The ordinates are in counts of each pixel. In Panels 1 and 3, the structures near 0 and 
650 km s“^ are night sky lines. In Panels 2 and 4, the structures at 600 km s“^ are also 
night sky lines. 

Figure 3a: The reduced data from one of our 17 subfields (the one just northwest of 
the SN). The range of data is 120 arcsec in the vertical direction (along the slit), with a 
free velocity range of 120 km s“^ corresponding to each of the hve slits in the horizontal 
direction. The SN itself is seen in [N H]A6583 in the lower (southern) left corner. 

Figure 3b: On the same scale as Figure 3a, the inferred velocity components are 
shown for the same data. The emission strength of each component (slit into 13.3 arcsec 
subsamples) is indicated by the width of each line sgment, which varies with the logarithm 
of linestrength. The tickmarks along the vertical axis indicate divisions of 13.3 arcsec and 
those along the horizontal indicate steps of 50 km s“^ (in the frame of a single slit). 
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Figure 4a: 3-D surface plot of all eight velocity components in front of SN1987A. They 
are collected into four groups, A, B, C and D (described in §4), with labels of each letter at 
the corners of each group indicating its approximate extent on the plot. The vertical axis is 
Vhei in km s“^. The other two axis are RA and Dec in J2000. North is to the left. West is 
the direction pointing into the paper. Note that structure within about 5 km s“^ will blend 
by one of them disappearing. 


Figure 4b: 3-D surface plot of the three velocity components in Group C, 265, 277 
and 285 km s“^. The axes and orientation are the same as in Figure 4a. Notice that the 
277 km s“^ component joins the 265 km s“^ component to the southeast of SN1987A, and 
the 285 km s“^ component to the south and slightly east of the SN. 


Figure 5: Flux contour of combined Vhei = 265, 277 and 285 km s“^ components 
on top of the R1170-complex image ( pCu, Grotts fc Kunkel l9^) . The axes are RA and 
Dec (J2000). Both the flux- and the dust-maps are smoothed in both directions over 13", 
which is the resolution limit of the echelle survey. This survey only covered the region 
within the dashed lines. The “holes” in the gray scale dust-map are due to the lack of data 
and inadequate contrast in image display after smoothing. The readers are encouraged to 
compare this plot and Figures 6 and 7 with the similar plot in |Xu, Grotts fc Kunkel 1995| . 

Figure 6: Flux contour of the Vhei = 235 km s“^ components plotted over the combined 
SE3140 and N3240 signals in greyscale (pCu, Grotts fc Kunkel 1995|) , analogously to Figure 

5. 


Figure 7: As in Figure 5 and 6, flux contour of combined Vhei = 245 and 255 km s ^ 


components plotted over greyscale image of the R430-complex (|Xu, Grotts fc Kunkel 1995|) . 
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